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ABSTRACT 

We will pursue an improved physical understanding and mathematical model for the boundary condition at an 
oscillating contact line at high Reynolds number begun by Ting & Perlin (1). We expect that the body force is locally 
unimportant for earth-based systems, and that the local behavior may dominate the mechanics of partially-filled 
reservoirs in the microgravity environment (Weislogel & Lichter 2). 

One important space-based application for this contact-line study is for Far ad ay- waves. Oscillations in the direction 
of gravity (or acceleration) can dominate the fluid motion during take-off and reentry with large steady-state acceler- 
ations and in orbit, where fluctuations on the order of 10“ 4 <7 occur about a zero mean. Our experience with Faraday 
waves has shown them to be “cleaner” than those produced by vertical or horizontal oscillation of walls. They are 
easier to model analytically or computationally, and they do not have strong vortex formation at the bottom of the 
plate. Hence many, if not most, of the experiments will be performed in this manner. The importance of contact lines 
in the microgravity environment is well established (Concus, et al. 3; Sigigner &: Weislogel 4; Capodanno 5; Kim, et 
al. 6; Lucassen, et al. 7; Merte, et al. 8). 

We will compare high resolution measurements of the velocity field (10//m resolution) using particle-tracking and 
particle-image velocimetry as the fluid/fluid interface is approached from the lower fluid. The spatial gradients in the 
deviation provide additional means to determine an improved boundary condition and a measure of the slip region, 
unavailable from the work by Ting h Perlin. Dissipation, the size of the eddy near the contact line, and hysteresis will 
be measured and compared to linear and nonlinear models of viscous and irrotational but dissipative models. 

1 INTRODUCTION AND THE STATE-OF-THE-ART 

The contact line is the intersection between two distinct fluids and a solid: here, the fluids are water and air, the 
solid is glass. The contact angle 9 C is defined as the angle between the tangent to the fluid interface at the contact 
line and the water-glass interface (see figure 1). Contact-line behavior may be couched as a relationship between the 
contact angle and the relative (liquid-to-solid) motion of the contact line. “Contact angle” is used in place of “apparent 
contact angle.” (the macroscopic contact angle, Dussan V. 9) To determine and quantify the contact-line behavior, 
precise measurements of the contact-line position (and thus velocity and acceleration) and dynamic contact angle (z.e. 
contact angle with a moving contact line) must be conducted. 

Information regarding the contact line and contact angle in uni-directional viscous flows abounds in the literature. 
A complete and thorough discussion of static contact angle and dynamic contact angle for uni-directional relative 
velocities is given by Dussan V. (9). These results are representative for experiments with small Reynolds number, 
small Weber number, and small Froude number. For these cases, there is no inertial effect from contact-line movement. 
The relationship shows that the contact angle is a function of the contact-line velocity only, i.e. 6 C — f(V r ). In figure 
2, a r represents the receding angle defined when V r — ► 0“ (the contact point moves toward the liquid), a a represents 
the advancing angle defined as V r — ► 0 + (the contact point moves away from the liquid). Because of the hysteresis, the 
contact-line boundary condition is nonlinear, even for low Reynolds number uni-directional, steady motion. 

A difficulty arises in the solution of the viscous flow field — even with a (uni-directionally) moving contact line. 
Dussan V. h Davis (10) showed that a non-integrable-stress singularity occurs at the contact line if a moving contact 
line is forced to obey the no-slip boundary condition. A popular model to avoid the singularity that allows the contact 
point to move along the solid surface was originated by Navier. He postulated that the resistant stress at a fixed solid 
surface is proportional to the relative velocity or ftu = fidujdn . Here, /? is a constant, u is the velocity parallel to 
the fixed solid surface, n is the distance normal to the solid surface, and /i is the liquid dynamic viscosity. A terse 
description of this slip boundary condition is given by Goldstein (11). 

As opposed to their uni-directional counterparts, oscillatory contact-line boundaries have not been studied exper- 
imentally in a significant way until very recently (1). Oscillatory contact-line boundaries are important in wave and 
other flow interactions with a solid boundary. Benjamin & Scott (12) and Graham-Eagle (13, 14) investigated waves 
propagating in a narrow open-channel with pinned-end (fixed contact-line position) edge conditions. Hocking (15) 
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demonstrated the importance of surface tension effects at a contact line and showed that the damping of gravity- 
capillary waves at a rigid boundary is mostly due to capillary effects. Several other publications that demonstrate the 
importance of capillary effects in wave motion and damping include Miles (16, 17, 18); Hocking k Mahdmina (19); 
Joo, Schultz k Messiter (20); Cocciaro, Faetti k Nobili (21); and Cocciaro, Faetti k Festa (22). 

Quantitative information is sparse regarding oscillating contact lines. Based on the relationship given by Dussan V. 
(9) for uni-directional contact lines, Young k Davis (23) proposed four possible relationships between contact angle and 
contact-line velocity for an oscillating contact-line boundary: (a) with contact-angle hysteresis; (b) fixed contact line; 
(c) fixed contact angle; and (d) smooth contact-angle variation (no contact-angle hysteresis). Using these relationships 
and a force balance at the contact line, they presented a solution to the oscillatory contact-line motion in the creeping- 
flow limit. They considered small motion of the plate so that the inertial effect is small; therefore, the contact-line 
behavior is governed by a relationship similar to that of Dussan V. (9). They found that contact-angle hysteresis and 
steepening of the contact angle with increasing contact-line speed are dissipative effects. Also, they stated that the 
contact-line motion tends to lag behind the plate motion due to inertia. 

Hocking (24) used an oscillatory contact-line boundary condition to calculate the waves generated by a vertically 
oscillating, vertical plate. Capillary effects were included in his analysis at the contact line. Two approximations 
originally proposed by Young k Davis (23) to the uni-directional model of Dussan V. (9) were used. Miles (17), 
addressing the same problem as Hocking, also used a Navier-slip boundary condition along the plate associated with 
a boundary condition at the contact line. In his analyses, viscosity was included and a non-zero initial free-surface 
meniscus was also considered in one case. A slip length, l 3 (defined as slip velocity divided by shear at the wall), was 
adopted to replace the no-slip boundary condition along the entire plate. By a phenomenological hypothesis, he posed 
the slip boundary condition (equivalent to Navier's slip model with l s = pi/ (3) 

V V -i dv 

V ~ V solid — , 

on 

where V is the vertical velocity along the plate (|| at the contact line). He stated that l s is a function of position 
along the plate surface and vanishes at a distance, l u — yj2v/u> (the viscous length scale), away from the contact line. 
That is, the flow obeys the no-slip boundary condition at some distance i n away from the contact line. For simplicity, 
he assumed that /, is constant, but experimental results show this is unrealistic as mentioned by Miles (17). 

From the discussion and conclusions of Ting k Perlin, the edge condition at an oscillatory solid boundary (i.e. a 
trijunction) is shown to be very important, but it is neither fully understood, nor predictable for large contact-line 
motion. Many theoreticians have adopted some form of the uni-directional model and a slip boundary condition along 
the solid surface, however, according to measurements by Ting k Perlin, this model is incorrect for large oscillatory 
(sinusoidal) plate motion. The qualitative behavior of the contact-line movement changes as the amplitude of the plate 
oscillation changes (figure 3, 4). The mean slip coefficient changes dramatically with oscillation amplitude as shown in 
figure 5. 

Dissipation due to viscosity (and the vortex imaged by Ting k Perlin) and hysteresis are very important, ubiquitous 
phenomena. We will accurately measure this flow field and compare it directly to numerical results. As an example 
of the lack of understanding of dissipation in oscillatory flows, many published results present a discussion stating 
that a one to three order-of-magnitude increase in the viscosity is required to match numerical results with measured 
dissipation (e.g. Hocking 15 and Henderson k Miles 25). We believe this so-called effective viscosity is the result of 
contact-line- generated dissipation . 

2 OBJECTIVES 

The main goal of this research is to provide a physical explanation of oscillatory contact-line behavior, a more accurate 
boundary model for large-amplitude oscillations, and an attempt to model dissipation in the vicinity of the contact 
line. The proposed problem clearly requires an energetic, multi-disciplinary approach. 

Specific objectives are: (1) to conduct precise measurements of the flow field in the immediate vicinity of the 
contact line and as a function of position along a vertically oscillating reservoir or upright plate using p article- image- 
velocimetry techniques; (2) to compare the measured velocity fields with those predicted by the exact analytic solution 
of Stokes, thus defining the slip region as well as elucidating the effect of the contact line, meniscus, and the interface 
on the local fluid physics, this by investigating the deviation from the exact solution; (3) to develop a proper boundary 
condition for the stick-slip motion of the contact-line region (Note that the Ting k Perlin model did not include the 
effects of acceleration directly.); (4) to develop a numerical code to simulate the flow field using the boundary condition 
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developed and the viscous equations of motion; and (5) to verify and quantify the accuracy of the boundary conditions 
developed by comparing the interfacial elevations and flow fields of the measurements with those of the numerics. In 
the following, we discuss the laboratory facility, experimental techniques, and the experiments; the numerical methods 
and simulation; the contact-line effects on Faraday waves; and the comparisons to be conducted. 

3 LABORATORY FACILITY, EXPERIMENTAL TECHNIQUES, AND EXPERIMENTS 

The study will include a synergistic set of physical and numerical experiments as follows: a contact-line model will 
be developed based upon physical experiments and the deviation of the experiments from the Stokes second problem 
solution; the newly determined boundary condition used in conjunction with the viscous equations and other boundary 
conditions will be solved numerically ; and a comparison of the physical and numerical wave fields generated will quantify 
the effectiveness and applicability of the contact-line boundary condition. In addition, the energy dissipation at the 
contact-line region will be studied, especially with regard to the relatively large vortex discovered by Ting k Perlin. 

Using a time series recorded by a high-speed camera of the contact-line region as illuminated by a laser sheet in the 
presence of dye-laden water, one obtains the position of the plate, the position of the contact line, and the apparent 
contact angle. Sufficient magnification is achieved with a proper choice of tenses and extension tubes; however, to obtain 
a detailed velocity field with a large number of velocity vectors still requires resolution available only on film. The 
relatively poor resolution of the video system was one of the limiting factors on the data of Ting k Perlin (1). These 
data had insufficient resolution to obtain meaningful accelerations through two differentiations of the plate/contact- 
line position. Our experience shows that 4”x5” Kodak Tmax 400 film works well for velocimetry with a film resolution 
on the order of 5 microns. Using the computer to generate the timing sequence, any desired phase of the oscillation 
can be captured on the film. In addition, a technique developed by Lin k Perlin removes light reflections present in the 
fluid-fluid interface, thus removing the ambiguity of the interfacial location and recording only those particles actually 
located in the lower fluid. The region of high shear (as seen by Ting k Perlin) immediately adjacent to the contact-line 
region will be studied also. From these data, we believe a reasonable boundary condition is achievable. We will have 
inferred the region (length) of slip through the comparison with the solution of the oscillating plate in a semi-infinite 
fluid. An improved boundary condition from that used by Ting k Perlin (1) will be formulated, since the acceleration 
information will now be available. 

Experiments to verify the boundary condition will be conducted to determine the interfacial elevation and the 
particle velocities some distance away from the plate. These data will be used for comparison to the numerical solution 
of the flow field obtained with the boundary condition developed herein. Thus, the validity and accuracy of the model 
will be determined. We also intend to evaluate and compare other independent measures of contact angle (Concus, et 
al. 3; Seebergh k Berg 26; Li, et al. 27; van der Zanden k Chesters 28). Many of these are adapted to the microgravity 
environment, but may not be applicable for oscillatory flows. 


4 NUMERICAL METHODS AND SIMULATION 


We will model the flow near the contact line in a 2 x 2 matrix of linear and nonlinear and inviscid and viscous flows. 
The inviscid model has been developed by Hocking (6), but only for linearized boundary conditions and with a simple 
contact-angle model that bears no resemblance to figure 2. On the other hand, there have been fully-viscous analyses 
by several groups (Kafka, et al. 29; Dussan, et al. 30), but for steady, creeping flow. 

The flow we wish to model is high Reynolds number. Hocking (15, 24) showed that many of the major features 
(such as damping) can be retrieved from an inviscid, even irrotational model. This is important as modeling very high 
Reynolds number flows with a fully-viscous computation is quite daunting. Features such as slip length and the eddy 
near the contact line must of course be discarded. 

We will extend Hocking’s analysis to include a more general model for contact-line behavior, namely 

*c = W,^). 


This model can be made sufficiently general to accommodate the experimental data of figure 3. This data shows 
that the contact angle versus relative speed has hysteresis itself, as opposed to the “usual” contact angle hysteresis 
represented by the vertical line in figure 3. A simplification when no contact angle hysteresis exists is possible for 
periodic motion with a complex A (Miles 17), 


V r = A 


drj_ 

dx' 
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at the expense of fitting the data of figure 3 only crudely. 

When we model the viscous flow, we again resort first to linearized flow to simplify the problem and make it more 
amenable to analysis. Now a slip model is required. One important aspect that is not usually studied is that the slip 
coefficient in the unsteady case must be compatible with the contact angle model to avoid reintroducing a singularity 
at the origin. Specifically, at the contact line, we have two conditions on slip — a contact-angle condition and a slip 
condition. For the simplest contact angle relationship we get 


V r 



_ . dVcL 
5 dx 


Since Vcl — dr]/dt , this requires that A = l s d/dt. More complicated compatibility constraints are required for more 
general contact angle relationships. 

In the last year or two of the contract period, we expect to run spectral element computations of modestly high 
Reynolds numbers. Then the linearized free surface boundary conditions can be relaxed and comparisons of eddy size 
and strength can be compared with the experiments. 


5 CONTACT-LINE EFFECTS ON FARADAY WAVES 

We review and discuss some pertinent results from Jiang, et al. (31) as regards contact-line effects on Faraday waves. 
Specifically, we discuss the competition between viscous effects (tend to decrease frequency) and contact-line effects 
(tend to increase frequency). The threshold forcing amplitude for the fundamental mode (wavenumber ki= 1) was 
measured and a stability diagram is presented in figure 6 for three slightly-different fluids and two tank-aspect ratios 
(AR). Figure 6(a) presents the measurements of water (hollow triangle, AR=10 and hollow diamond, AR=2.61) and 
mixed water with fluorescein dye (hollow square, AR=10). Also shown is the subharmonic Faraday wave with two 
wavelengths in the tank (solid triangle, 2, AR=10). Figure 6(b) shows the results for water mixed with Photo Flo 
200, volume ratio 100:1 (solid square, AR—2.61 and solid diamond, AR=10). Dimensionless frequency p (= 
where u>i is the frequency of the excited mode and is the forcing frequency) and dimensionless forcing amplitude 
q (= 2fki tanh k{h where f is the forcing amplitude and h is the water depth) are evaluated based on the inviscid 
natural frequency If p and q lie within the upper and lower solid diagonal lines (unstable region), then subharmonic 
resonance will occur according to linear theory (Benjamin & Ursell 32); otherwise, the liquid surface is not excited 
(stable region). The dashed curves represent weakly-nonlinear predictions of Henderson & Miles (24), now denoted by 
H&M. The hollow circles represent their experiments (with Photo Flo present), shown in figure 6(b). 

According to H&M, the viscous natural frequency u> v should be less than the inviscid natural frequency u,-, i.e. 
Pmin j the p value associated with the minimum of q (g mjn ), should be larger than one. However, our experimental 
results show an increase in u> v , not a decrease. As shown in figure 6(a), the p m i n associated with each set of symbols is 
less than one, implying that u v (= | 9min ) > u>i . Benjamin & Scott (13) measured waves propagating in a brimful 
channel and concluded that the contact-line effect increases the wave celerity (phase speed). By reducing the tank 
aspect ratio and lessening the relative effect of the contact line by the addition of Photo Flo, we demonstrate this same 
phenomenon is responsible for our increased u> v for fixed wavenumber. 

The aspect-ratio effect on the frequency is seen in figure 6(a) by comparing measurements in the wider tank 
(AR=2.61) and those in the narrow tank (AR=10). Due to weaker contact-line effect in the wider tank, u) v should be 
and is smaller, corresponding to an increase in p m in ■ Also, since the contact line and wall-boundary layers contribute 
the most to the dissipation, the effective damping in the wider tank is less owing to the smaller surface-area- to- volume 
ratio. Therefore </ min in the wider tank is smaller. 

Figure 6(b) represents our measurements with the mixture of water and Photo Flo. Since the addition of Photo Flo 
reduces the contact-line effect, p mm and q m i n are relatively more affected by viscous damping and the neutral-stability 
curve shifts towards HfcM’s prediction. H&M’s theory still overpredicts q on the lower side of the neutral-stability 
curve and underpredicts u>„. Their experiments with a water-Photo Flo mixture (hollow symbols) agree with their 
theory quite well. 

6 COMPARISONS 

Measurements and numerics will be compared in detail. We will have a “closed-loop” prediction in the sense that 
our experiments will determine boundary condition along the plate. We will generate numerical solutions with these 
boundary conditions and finally compare them to the flow field adjacent to the plate to quantify and verify the 
applicability of our boundary condition. 
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Figure 1. Sketch of contact-line region. 
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Figure 2. Contact angle vs. relative velocity. 
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Figure 3. Typical graph of the contact angle vs. Figure 4. A pictorial sketch of particle trajectories 
nondiinensional relative velocity. near the contact line. 



Figure 5. Mean slip length versus stroke 
amplitude. 



Figure 6. Subharmonic instability limit curves of 
Faraday waves for several experiments. 
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